Critical wave functions in disordered graphene 
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In order to elucidate the presence of non-localized states in doped graphene, an scaling analysis 
of the wave function moments known as inverse participation ratios is performed. The model 
used is a tight-binding hamiltonian considering nearest and next-nearest neighbors with random 
substitutional impurities. Our findings indicate the presence of non-normalizable wave functions 
that follow a critical (power-law) decay, which are between a metallic and insulating behavior. The 
power-law exponent distribution is robust against the inclusion of next-nearest neighbors and on 
growing the system size. 
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Graphene is a two-dimensional atomic crystal 1 1 
with the highest known charge carrier mobility 2| 
and thermal conductivity Q at room temperature. 
Both properties set up graphene as a raw material 
for transistors, however, the ' graphenium inside' era 
is quite far [J]. Keeping in mind the aim of de- 
sign a transistor, the problem turns out to be how 
to alchemize it into a semiconductor. An alterna- 
tive is to dope graphene. This leads immediately 
to the question of having quantum percolation in 
two dimensions, which has been the subject of de- 
bate many years ago Q . In the literature, usually 
it is found that "there is no true metallic behavior 
in two dimensions" as a consequence of the fact that 
all eigenstates are localized even when the disorder is 
weak . I n grap hene, there has been a debate about 
this point [8l4lOj . Recently, disordered graphene has 
been classified using symmetry arguments around 
the Dirac point [ll[ ; this classification allows a min- 
imal conductivity behavior. Experimentally, the 
minimal-conductivity was measured for potassium 
atoms onto the graphene surface 12, 13|. Further- 
more, it has been found that a metal-transition 
can be observed when graphene is doped with H 
14j . Also, using a non-interacting electron model 
enriched by first-principles calculation the metal- 
insulator transition is observed [ljl [lj|. In this 
Letter, we present numerical evidence that shows 
a very interesting scenario. We have characterized 
the probability distribution of the moments associ- 
ated to the wave function using the inverse partic- 
ipation ratios. The scaling of this quantity is fre- 
quently used to discriminated when an cigenstate 
is extended or localized. In doped graphene, we 
found extended states which do not follow the usual 
exponential localization; instead, these are critical, 



i.e., the wave function decay spatially as a non- 
normalizable power-law. This behavior evidenced 
the multifractality of the wave function [T3] ■ Notice 
that in the orignal development of the scaling theory, 
critical states were not considered Q ■ 

As a model we use the tight-binding hamiltonian, 
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where the nearest neighbor, t = 2.79 eV, and next- 
nearest neighbor (NNN), t' = 0.68 cV, hoping pa- 
rameters are included; these values have been taken 
from reference [l8j ]. The impurity sites, £, have been 
distributed randomly in the lattice according with a 
fixed concentration, C; and e is the self-energy for 
an impurity. 

In order to investigate localization, we introduce 
the inverse p-participation ratios (IPRs), 
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where \&k is the wave function associated to the 
eigenstate k with energy E^, which solves the 
Schrodingcr equation H^u = E^^. The index i 
belongs to the sum over sites, N is the total num- 
ber of sites and p is an integer. When p = 1, 
||^k||2 = 1 because of the normalization condition. 
If the wave function of the eigenstate follows the 
power-law, ^E , (r)k ~ |r| _Q , the p-IPRs are scaled 
as [H, 
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when p > 1. The behavior TV -1 corresponds to a 
metal while N° corresponds to an insulator. Notice 
that here N cc L 2 . where L it is the length sample. 

To evaluate Eq. (J3]) , we calculated all the eigenval- 
ues and eigenvectors of TL by numerical diagonaliza- 
tion. The p = 2 IPR behavior is shown in Figurc[T]as 
a function of N for several selected energies, using 
different impurities self energies without and with 
NNN interaction. For the energy E c = E^ — OAt, far 
from the Dirac energy, we compare the behavior for 
pure graphene [(green online) diamonds] with doped 
graphene [(blue online) squares] . For pure graphene, 
the state is extended since the p=2 IPR goes like 
TV -1 , as shown in Figure [TJ For doped graphene, 
is clear that the p=2 IPR can be fitted with a line, 
which suggests a non-localized, power law behav- 
ior. On the other hand, for states near the Dirac 
energy, in doped graphene the (dark blue online) 
circles and (red online) triangles behavior suggests 
localized states, since the p=2 IPR scales as oc TV . 
Thus, from this scaling analysis is clear that doped 
graphene, even in the absence of NNN interaction, 
presents different kinds of localization, as has been 
suggested in reference [2(| due to frustration effects 
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as well as in experiments [14| . 
In order to know the exponent distribution of the 
power law behavior, we introduce the integrated dis- 
tribution of exponents fl9| , 



4 8 



(4) 



where O is the step function. The exponent distribu- 
tion is plotted in Figure [2] for several impurity types 
and two different concentrations, 1% and 5%. For 
pure graphene, all the states have the same scaling 
behavior, and a step is observed at 7 = —3. This 
means that all states have the same scaling, and thus 
all are extended, as expected from Bloch's theorem. 
However, for doped graphene, we observe two main 
effects. First there is a shift to higher values of 7 
and second, the jump is not anymore a discontinu- 
ity. Instead, we observe states that have a distri- 
bution of 7 values. In all cases, we observed that 
the minimal value of 7 is ~ —5/2, which means that 
the most extended states, follow a power law that 
goes as r -3 / 8 . States with 7 m are exponentially 
localized. Since no clear jump is observed in the val- 
ues of 7, it seems that there is a range of values for 
the exponents of the critical wave functions. Also, 
it is worthwhile mentioning that the presence of the 
NNN interaction preserves this behavior, which al- 
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FIG. 1. (Color online) The p = 2 IPR behavior as a 
function of N for several selected energies, using different 
impurities self energies (from top to bottom) and with- 
out and with NNN interaction (left and right columns, 
respectively). The selected energies are the Dirac point 
energy [En, (dark blue online) circles], the resonant state 
energy [E T , (red online) triangles] and E c = Er> — OAt 
[(blue online) squares] for doped graphene ' d '. In all 
cases, doped graphene has a 5% impurity concentration. 
For the state with energy E e , we present in all cases 
the scaling exponent that results from the fitting, shown 
in the figures with (blue online) lines. For comparison 
proposes, we include the case of pure graphene for the 
energy Et [(green online) diamonds]. 



lows its experimental verification since the NNN is 
always present. 

To verify that such behavior is preserved as the 
system grows, in Fig. [3] we present the distribution 
1(7) for different sample sizes. We can observe that 
the behavior is very similar at all sizes. From this 
analysis, we can conclude that there are many non- 
cxponcntially localized states, and that the power 
law behavior is not a finite size lattice effect. From 
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FIG. 3. (Color online) Example of the Integrated dis- 
tribution of exponents for different sample sizes using a 
fixed concentration of impurities (5%) with self-energy 
e — —6t. The left panel corresponds to the model with- 
out NNN interaction, and the right panel graph with 
NNN. 



FIG. 2. (Color online) Integrated distribution of ex- 
ponents for different impurity configurations using lat- 
tices with TV = 16, 928 sites. Extended states are at 
7 = — 3 and localized at 7 = 0. We observe that for 
pure graphene [(green online) diamonds], the distribu- 
tion is a step function at 7 = —3, meanwhile for doped 
Jiapheiie Lire disQibuliun is shifted and H is iiu lunger a 
i top function. 



the values of 7, these atatca arc critical and non 



hormalizablc. Finally, a careful check of auch atatca 
reveals that localized states are near the Dirac point 
and at the band edges, while the power law non- 
normalizablc states are near the middle part of the 



valence and conduction band, in agreement with prc- 
vious theoretical arguments 



In conclusion, using a scaling analysis of the par- 



ticipation ratio, we have shown that the presence of 
disorder in graphene do not localize exponentially all 
states; instead, some states are critical with a dis- 
tribution of exponents. This result is robust against 
the inclusion of NNN interactions, in which the chi- 
rality is not preserved. This result is not only im- 
portant for graphene, but it makes a revival of an 
old discussion concerning the possibility of having 
lanumaluus Liuaiiliun pcieulaLiun in two dimensional 
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